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The two-orbital model for manganites with both non-cooperative and cooperative Jahn- Teller 
phonons is studied at hole density x=0.5 using Monte Carlo techniques. The phase diagram is 
obtained varying the electron-phonon coupling and the t2 g -spins exchange. The insulating CE-type 
charge- and orbital-ordered state with the z-axis charge stacking observed in narrow-bandwidth 
manganites is stabilized in the simulations. Its charge gap Aco is much larger than the critical 
temperature kaTco- Metallic-like A-type and ferromagnetic states are also obtained in the same 
framework, and the phase boundaries among them have first-order characteristics. 



The complicated interplay between charge, spin, or- 
bital, and lattice degrees of freedom is believed to in- 
duce the unexpected magnetic and transport phenom- 
ena observed in Mn-oxides, such as the 'colossal' mag- 
netoresistance (MR) effect ||]. The subtle properties of 
manganites are especially complex at 50% hole density 
where experiments have established that an insulating 
charge-ordered (CO) state exists in Ndo.sSro.sMnOa 
Pro.5Cao. 5 Mn0 3 g, LaSr 2 Mn 2 7 @,and several oth- 
ers. The state has alternating d3 X 2_ r 2/d3 y 2_,.2 orbital 
and CE-type antiferromagnetic (AF) order Its melt- 
ing by a magnetic field leads to a huge negative MR || , 
probably induced by mixed-phase tendencies 

In spite of its importance, the origin of the x=0.5 CO- 
state is still unclear. The observation of charge popu- 
lating only the even (or odd) Mn-ion sublattice in the 
xy plane naively suggests that the nearest-neighbor (NN) 
Coulomb repulsion Vnn is important for its stabilization. 
However, the z-axis stacking of charge |9| and existence 
of bistripes at x>0.5 [(To), both penalized by a strong 
NN-repulsion, show that Vnn is smaller than expected 
and other ingredients are needed to understand the CO- 
state jnj. In addition, at x=0.5 experiments exhibited 
two other phases in competition: (i) A d x 2_ y 2 orbital- 
ordered (00) but charge-disordered (CD) state with A- 
type AF state and two-dimensional (2D) mctallicity (in, 
e.g., Pr .5Sr .5MnO3 g]). (ii) A metallic spin ferromag- 
netic (FM) CD and orbital-disordered (OD) state (in, 
e.g., Lao.5Sr .5Mn03). For a dominant Vnn, stable CD 
states at x=0.5 are difficult to understand. It may oc- 
cur that electrons coupled only through on-site Coulomb 
interactions produce this exotic behavior, as suggested 
by Hartree-Fock calculations [[l2|. However, in this con- 
text unbiased many-body calculations are difficult and 
computational studies have technical complications such 
as the sign-problem. In addition, the 'colossal' oxygen 
isotope shift of the CO transition in Ndo.sSro.sMnOs 
(NSMO) jl3| cannot be explained using purely Coulom- 
bic approaches, which predict no isotope effect. 

Searching for an alternative mechanism to stabilize 
the charge-stacked CO-state of manganites, in this pa- 
per the first comprehensive computational analysis of the 



half-doped two-orbital model with both non-cooperative 
and cooperative Jahn- Teller (JT) phonons in the large 
electron-phonon coupling (A) regime is presented. Our 
main result is that a CE-type CO orbital-ordered ground 
state can be stabilized in this framework with charge 
properly stacked along the z-axis. Many other features 
found in the simulations are also in excellent correspon- 
dence with experimental observations. 
The two-orbital Hamiltonian is given by 

H =- X] ^bcL^i+vbo— Jh^S; • Si + JAF^S; • Sj 

ivtrab i (ij) 

+ X Y1 C L(Q2i^l + Q3iCT 3 )abC ibff + ^2(Qli + Qli), (i) 
iabcr i 

where cn a (ci2 CT ) is the destruction operator for an e g - 
electron with spin a in the d x 2_ y 2 (d 3z 2_ r 2) orbital at 
site i, v is the vector connecting NN sites, and t^ b is 
given by t^ c 1 =— V3t^ 2 =— v / 3t 2 i 1 =3t2 < 2=t (energy unit) for 
v=x, t^^v^t^ =v / 3t 21 =3t 2 ' 2 =t for v=y, and t| 2 =4t/3, 
tf 1 =tf 2 =t| 1 =0 for v=z. The Hund coupling J H (>0) 

links the e g -electron spin SF=X)aa/3 c iaa <Ta /3 Cia / 3 anc ^ tnc 
calized t2 g -spin Si assumed classical (] Si |=1) . The eg- 
electron density (n) varies with a chemical potential /i 
(the hole density is x=I— (n)). Jaf is the AF coupling be- 
tween NN t2 g -spins. The fourth term couples e g -electrons 
to JT-mode distortions, Q2i and Q 3 i, assumed classical 
HQ- o~ i and cr 3 are Pauli matrices. When the {Q} are 
independent, "non-cooperative" phonons are used |p|,pd|. 
In real materials, however, adjacent MnOe octahedra 
share an oxygen. In this "cooperative" case, the proper 
variables are the oxygen displacements u^ from equilib- 
rium along the Mn-Mn bond in the v-direction, and ex- 
tra terms are added to Eq. (I) for the breathing- mode 
phonon [l5| . Standard Monte Carlo (MC) simulations 
for the classical spins and phonons are here used 

Figure 1(a) already illustrates one of the main re- 
sults of the present study. Shown are the MC real- 
space charge correlations at x=0.5 for representative cou- 
plings in 2D at low temperature(T). For intermediate Jaf 
(0.1 < Jaf 0.2), the charge correlations were found to 



1 



be robust at all distances and they are positive (negative) 
on the even (odd) sites, compatible with the CO-pattern 
observed in x=0.5 experiments. Our results show that 
this behavior is essentially independent of A for interme- 
diate Jaf, although the CO-pattern is more easily ob- 
served for large A. For small Jaf (^S 0.1), on the other 
hand, the CO-state appears only for large A. The ob- 
served CO-state originates from the tendency of the JT- 
distorted Mn 3+ -ions to maximize their relative distances 
to exploit the kinetic energy of the mobile carriers. Fig- 
ure 1(b) contains the Fourier transform of the t2 g spin 
correlations S(q) in the CO-regime. For intermediate 
Jaf, peaks at q = (n, 0) and (7t/2,7t/2) are observed in 
excellent agreement with the result of a perfect CE-spin 
arrangement. Note that size effects appear small in these 
calculations, and throughout the paper. For small Jaf, 
only one peak is detected at q = (0,0), indicating the 
FM phase. Under the CE-type spin environment, it can 
be shown that the CO-state is stable even for small A, 
since the system is effectively one-dimensional (ID) due 
to charge confinement along the zigzag FM-chain lfll. In 
this sense, Jaf plays a role more important than A for 
the CE-state stabilization. This is in contrast to the FM 
state where a finite and large A is needed to stabilize the 
CO-state due to its higher dimensionality 



and 2.0. For small Jaf, #a = &b = 7r(d x 2_ y 2), since the 
orbitals prefer to have large overlaps both in the x- and 
y-directions. On the other hand, for intermediate Jaf, 
9 a = 27r/3(d 3x 2_ r 2) and 9b = 47r/3(d3 y 2_ r 2). To under- 
stand this result, note that the orbital is polarized along 
the hopping direction to improve the kinetic energy. Due 
to this effect, the zigzag hopping path in the CE-type CO- 
state induces the alternatingd3 X 2_ r 2/d3 y 2_ r 2 orbitals. The 
above results in charge, spin, and orbital correlations for 
intermediate Jaf are schematically shown in Fig. 1(d), 
which is the CE-type CO-state observed in experiments. 
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FIG. 1. (a) Charge correlation C(m) = {nini4- m )-(n) 2 
with nj=y^ c iao- CiaCT f° r Jh=oo, T=l/100, and x=0.5, us- 
ing a 4x4 cluster with non-cooperative phonons and periodic 
boundary conditions. The squares (circles) [triangles] were 
obtained at J AF =0.1, A=2.0 (0.05, 2.0) [0.05, 1.6]. (b) Spin 
structure-factor S(q)=J^ e lq ' m (Si ■ Si+ m ) for the t2 g -spins 
at the same parameters of the squares in (a). The squares 
(circles) here denote the 4x4 MC (ideal CE-type state) re- 
sults. Triangles are MC results on an 8x8 cluster, (c) 8 a and 
0b (see text) vs. Jaf for A=1.6 and 2.0. (d) Charge, spin and 
orbital arrangement compatible with our MC results. 

Regarding the orbital arrangement, let us consider 
sites with substantial charge (A and B in Fig. 1(d)) 
|^7| . The occupied orbital at each site is given by |2)i = 
-sin(fV2)|l)+cos(fV2)|2) with 9i = tan- 1 ((Q 2i y(Q 3i )). 
In Fig. 1(c), 9 a and 9b vs. Jaf are plotted for A=1.6 



FIG. 2. (a) MC energy per site E/N vs. Jap at 
x=0.5, A=1.5, T=l/100, and Jh=oo with non-cooperative 
phonons. The FM state is charge-disordered and or- 
bital(uniform)-ordered. "AF(2)" is explained in Ref. Jl8[ . 
"AF" is AF in both directions. The straight lines correspond 
to the finite Jaf correction to Jap=0 MC results, when a per- 
fect spin configuration is assumed, (b) Drude weight Dw vs. 
Jaf at T=l/100, x=0.5, using a 4x4 cluster. The results with 
cooperative (non-cooperative) JT-phonons were obtained at 
A=1.5 (1.2). (c) Phase diagram of the 2D two-orbital model 
with non-cooperative phonons at low-T. The phase bound- 
aries were obtained as in (a) using 4x4 and 8x8 clusters. All 
transitions are of first-order. In the FM-phases, the OO-state 
is uniform in the d x 2_ y 2 orbitals, while in the CE-phase the 
pattern Fig. 1(d) was observed. 

Figure 2(a) shows a typical MC energy vs. Jaf for 
A = 1.5. At small Jaf, a FM-state is stabilized since even 
at large A the optimal background for electrons has all 
spins aligned. However, its energy is penalized by the AF 
exchange and at realistic values of Jaf an equal AF-FM 
mixture such as the CO CE-state minimizes the energy. 
The discontinuity in d(E/N)/dJAF signals a first-order 
transition, as in experiments. With increasing Jaf, the 
CE-state is destabilized by a novel state labeled "AF(2)" 
pSj . The metallic vs. insulating character of the vari- 
ous phases can be analyzed with the Drude weight Dw 
f§. Results at A=1.2 (Fig. 2(b)) show that the FM- 
state is metallic (Dw 7^0) and the CE-state is insulat- 
ing (D\y=0). The jump in Dw again signals a first-order 
transition. Based on results similar to those in Figs. 2(a) 
and (b), the phase diagram in 2D and low-T was con- 
structed (Fig. 2(c)). It presents a rich structure includ- 
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ing FM and CE-phases. In particular, it contains a CD 
FM-state with uniform orbital order, which after consid- 
eration of three dimensional (3D) clusters will become a 
fully A-type state with charge and orbital characteristics 
as observed experimentally S. Then, all the currently 
known states at x=0.5 exist in our phase diagram, and 
they are in close competition. It is conceivable that small 
temperature, doping, or lattice spacing changes may alter 
their balance, as found experimentally. It is also possi- 
ble that the other novel phases in Fig. 2(c), especially the 
CO and 00 FM-state, may appear in future experiments 
under suitable conditions. 
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FIG. 3. (a) Spin structure factors S(0,0) and S(7r,0) and 
charge structure factor N(7r,7r) vs. T for a 4 x 4 cluster with 
A=1.2 and Jaf=0.15. The high-T background, corresponding 
to on-site correlations only, is subtracted, (b) Same as (a) 
but at A=1.9 and Jaf=0.10. (c) DOS vs. u for 4 x 4 and 
8x8 clusters for Jaf=0.15 and X—1.2. (d) Schematic phase 
diagram in the (Jaf, T) plane for the small A (left) and large 
A (right) regions. PM means paramagnetic phase. 

Here the cooperative effect on the above results is 
briefly discussed. When the oxygen positions {u^}, not 
{Q}, are used as variables in the MC simulations, the 
charge, spin, and orbital correlations are found to be 
similar to the non-cooperative results. In particular, a 
clear evidence of a CO CE-type state was observed in 
a wide range of the phase diagram and to stabilize this 
state Jaf was found to play an essential role as in the 
non-cooperative case. For Jaf=0, a CO-state is still ob- 
tained at large A, but it is spin FM. In addition, note that 
the A-type phase of Fig. 2(c) has Q3 modes with an uni- 
form mean-value at every site. This amount to a static 
deformation of the oxygen octahedra, shortening the c- 
direction (or elongating the a- and b-directions), similarly 
as in experiments pp|| . If appropriate lattice constants 
are not used for cooperative phonons, the A-type phase 
may not be stabilized. For this reason studying both 
non-cooperative and cooperative phonons is important. 

It is instructive to calculate the T-dependence of the 
various order parameters to compare with experiments. 
In Fig. 3(a), S(q) at q=(0,0) (FM) and (tt,0) (CE) and 



charge structure factor N(7r, it) vs. T are shown for 
an intermediate A. At T~0.1t, S(0,0) grows with re- 
ducing T as if FM were to become stable. However, 
a t Tco^TN~0.05t a sudden transition to a CO CE- 
state is found (T CO ~250K for t~0.5eV). Note that the 
appearance of ferromagnetism in a temperature win- 
dow above Tco is in quantitative agreement with re- 
sults observed in Lao.sCao.sMnOs (LCMO) [ fl9[, and in 
(Lai_ z Nd z ) .5Sro.5Mn0 3 (LNSMO) for z>0.6@. The 
results at large A are qualitatively different (Fig. 3(b)). 
Here Tco >Tn, and only a weak FM signal is observed. 
These results are in good agreement with experiments 
for Lao. 5 Sri. 5 Mn0 4 @ and Pro.5Ca . 5 Mn0 3 @- No 
metallic FM-state exists in these compounds, and our 
analysis predict that they have a larger A than LCMO 
and LNSMO The density of states (DOS) at sev- 

eral T's are in Fig. 3(c). A robust charge gap Aco^0.9t 
at T~0.01t is clear, as well as a precursor pseudogap 
at higher T's. The ratio Aco/kBTco^l8 is remarkably 
close to the tunneling results (~19) for NSMO ||. 

Let us understand intuitively the several transition 
temperatures. In the small A region, the FM transition 
temperature Tc scales with t due to the double exchange 
mechanism, while Tn is determined by Jaf irrespective 
of A. In the FM phase, a CO-state never occurs in this 
region. However, if a sharp transition to the CE-type 
AF state occurs at T=Tn, the CO-state is also gener- 
ated, since the CE-type state is a ID band-insulator (ll| 
and it can be shown using a standard mean-field approx- 
imation that kBTco~£Fe _1 ^ A , where e-p is the Fermi en- 
ergy much larger than Jaf- Thus, Tco ~ Tn holds for 
small A. Note that the ratio Aco/kEjTco is large be- 
cause Aco^t in the band- insulating region. When A is 
increased, Aco and Tco smoothly increase, keeping the 
relation Aco^keTco- Eventually, at large A Tco over- 
comes Tn- In this regime, Tc is determined by the su- 
perexchange mechanism as Tc~t/A 2 . The above estima- 
tions are summarized in Fig. 3(d). In the small A region, 
the metallic-like FM phase appears irrespective of Jaf at 
high T, but at low T, the CE-type AF state occurs for 
intermediate value of Jaf and the CO-state is concomi- 
tant only to the CE-phase. In the large A regime, the 
CO-state occurs in the whole region for T<Tco> higher 
than Tc and Tn- At low T, the system becomes FM- 
insulator or CE-type insulator depending on Jaf- All 
these results are supported by the MC simulations. Note 
that previously it was widely believed that modifications 
in the tolerance factor tF by chemical substitution lead 
to bandwidth changes sufficiently strong to drastically 
affect the properties of manganitc compounds. However, 
our results established that the tF-induced changes in 
Jaf are equally important, and they should be carefully 
considered when studying Mn-oxides. 

A difficult to understand feature of the manganite CO- 
state is the z-axis charge stacking. To address the pre- 
diction of JT-phonon studies, in Fig. 4(a) the MC en- 
ergy of a bilayer cluster vs. Jaf with non-cooperative 
phonons is shown. At small Jaf the FM-state dominates, 
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as in 2D (Fig. 2(c)). The CE-state in its "Wigner Crys- 
tal" (WC) and "Charge Stacked" (CS) versions (with a 
N(q) peak at n) and (n,n,0), respectively) appear 

as excited states at Jaf=0. WC has lower energy since 
the electronic kinetic energy improves when the Mill- 
ions spread apart. However, with increasing Jaf the sit- 
uation rapidly changes since the CS charge arrangement 
(CE-planes, one over the other, with opposite spins) has 
all z-axis bonds properly antiferromagnetically aligned. 
On the other hand, in the WC arrangement, with one 
plane shifted from the other by a lattice spacing along x 
or y, half the z-axis bonds are AF and the other half FM, 
making its energy JAF-hidependent. For this reason, at 
realistic Jaf, the CS-state dominates. Since similar re- 
sults are obtained with cooperative phonons and cubic 
clusters (Fig.4(b)), the effect appears robust and inde- 
pendent of fine details. The important interaction to 
stabilize the CS-state is the z-axis Jaf exchange. 




FIG. 4. (a) Energy per site vs. Jap at A=1.5, Jh=oo, 
and T= 1/100 with both cooperative and non-cooperative 
phonons. Solid (open) circles are MC results using a 4x4x2 
(8x8x2) bilayer. Squares at Jaf=0 are the MC results on the 
4x4x2 cluster assuming the t2 g -spin background according to 
the convention in (b), in a standard notation. The straight 
lines were obtained as in Fig. 2(a). (b) Same as (a) but for a 
cube 4 3 . Note that with non-cooperative phonons the A- type 
phase is stable in bilayers and cubes in a narrow JAp-window. 
(c) Mean-field energy per site vs. Vnn for a bilayer. Cou- 
plings are realistic (U' is the on-site inter-orbital repulsion). 
Since Jh is assumed infinite, the on-site intra-orbital repulsion 
is irrelevant. Solid and dashed curves denote the CE-type CS- 
and and WC-structures, respectively. 

What is the influence of the Vnn repulsion in our re- 
sults? It will certainly penalize the CS structure, and 
the WC-state will eventually become the ground state as 
this repulsion grows. However, the CS-state has a better 
magnetic energy than the WC-state. Comparing the en- 
ergy per site gained (-Jaf) and lost (Vnn /2) in the 3D CS 
structure, a critical coupling V£j N =2JAF~0.15eV is esti- 
mated if JAF~0.15t and t~0.5eV. A more sophisticated 
mean-field approximation |p4f at realistic manganite cou- 



plings provides V^ N -0.22eV (Fig. 4(c)). Although the 
Vnn bare value is ~3.6eV, it is reduced by the dielectric 
constant e w 20-45 leading to V NN ~0.08-0.18eV for 
manganites. Thus, it is concluded that the manganite 
NN repulsion is weak enough to allow for the CS-state to 
be stable due to the effect of Jaf- 

Summarizing, evidence was provided that the CE-type 
CO-state with Aco /keTco^l and charge stacked along 
the z-axis observed at x=0.5 in several manganites can 
be stabilized using JT phonons. The competing FM and 
A-type states also appear in the simulations. Our results 
have established that purely Coulombic approaches are 
not the only procedure to stabilize CO-states in man- 
ganites, but strongly coupled electron JT-phonon sys- 
tems provide an alternative framework where controlled 
many-body calculations are possible. 
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